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RESONANCE TYPE INSTABILITIES 
IN THE GASEOUS DISKS 
OF THE FLAT GALAXIES 
I. The Acoustical Resonance Type Instability and 
the Absence of Vortex Sheet Stabilization on Shallow Water. 

J.V. Mustsevaya and V.V. Mustsevoy 

Linear analysis of vortex sheet stability in the rotating gaseous disk or shallow water layer 
shows that presence of a central reflecting surface changes system stability significantly. An effect 
of absence of vortex sheet stabilization has been fonnd as compressibility exceeds Landau criterion. 
The properties of mnltimode short-scale instability of aconstical resonance type are investigated 
and probability of its influence upon experiments on the rotating shallow water is discussed. 

Introduction. 

One of urgent problems of astrophysics is clearing of mechanisms, resulting occurrence 
of a wide spectrum of instabilities in differentially rotating gaseous disks. Given question 
is most heavily studied in connection with attempts of an explanation of an angular mo¬ 
ment tap from accretion and protoplanet disks, as well as reasons of generation of spiral 
structure of flat galaxies. The important place in these researches occupies consideration 
of hydrodynamic instabilities, caused purely by differentiality of rotation, that is velocity 
kink. For non-dissipative supersonic axisymmetic currents of non-selfgravitating gas three 
mechanisms are known supporting these instabilities for disturbances with wave vectors 
without 2 ;-componenet in cylindrical coordinate system, namely centrifugal one (Moro¬ 
zov 1977, Morozov 1979, Fridman 1990, Nezlin & Snezhkin 1990), an interaction of modes 
with opposite signs of energy density and resonant waves amplification on corotation radius 
(Papaloizou & Pringle 1985, Glatzel 1987, Savonije & Heemskerk 1990). 

At the same time, as the real astrophysical systems are characterized by significant 
radial density gradients, without any reasons one more version of the resonant mechanism 
of waves amplification is ignored — because of repeated reflection from area of the specified 
gradient and from a corotation circle with amplification on amplitude and constructive 
interferention of reflected waves (i.e. actually we are dealing unstable waveguide modes). 
For plainparallel currents this mechanism was analyzed in many works (see, e.g., Blumen, 
Drazin & Billings 1975, Drazin & Davey 1977, Morozov & Mustsevoy 1991, Morozov, 
Mustsevaya & Mustsevoj 1991). 

In the present work using extremely simplified model it is shown that the latter am¬ 
plification mechanism is effective enough — so that the instability caused by it is not 
suppressed even by centrifugal stabilization of disturbances. Thus in considered model the 
role of an internal reflecting surface plays not a sharp gradient of density, but firm wall. On 
the one hand, this significantly simplifies the dispersion equation and its analysis, on the 
other hand, the results received in such model can have the direct application to analogue 
laboratory modelling of spiral structure of galaxies generation on installation with rotating 
shallow water (see Nezlin & Snezhkin 1990 and quoted there literature). In a central part 
of this installation a cylindrical rigid surface (protective casing, interfering penetration 
of a working liquid in the bearing), is located; we discuss an opportunity of influence of 
disturbances reflection upon stability of the vortex sheet on shallow water. The situation 
with reflection from density jump will be considered in the subsequent parts of the article. 

1. Model and basic equations 
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It is investigated a dynamics of fonrier-harmonics of small distnrbances in snch form: 


f = f{r, z) exp(im(p + iut), 


where uj is complex freqnency, m = 0 , 1 , 2 ,... is nnmber of mode on an azimnth, imposed 
on eqnilibrnm system parameters, representing stationary differentially rotating disk-like 
confignration of polytropic non-dissipative non-selfgravitating gas with eqnation of state 
Pq = Poc'll'y {cs is adiabatic sonnd speed, 7 is adiabatic index). The model stability 
is provided by radial balance of centrifngal force and gradients of pressnre and external 
gravitational potential Tq and hydrostatic balance along ^-coordinate. 

We consider the solntions of eqnations averaging along 2 ;-coordinate (on opportnnity 
of snch approach see Part II). The eqnations describing dynamics of small distnrbances 
in snch system are eqnivalent to those in rotating shallow water accnrate to replacement 
p on liqnid layer depth i7, Cg on ^/gH and for 7 = 2 (Fridman 1990, Landan & Lifshits 
1988). Thns the strnctnre of installation bottom respects a radial strnctnre of gravitational 
potential. 

The angnlar rotation speed gets ont as 0(r) = Vtin at r < 0(r) = Vtex at r > 

and hrm reflecting snrface sitnated at r = Rp (eqnivalence of a hrm wall and density 
jnmp with signihcant gradient is shown by Morozov & Mnstsevoy 1991). It is necessary 
to note especially that the discontinnons model was chosen qnite conscionsly, in order 
to exclnde from consideration all instability mechanisms excepting considered. Really, for 
resonant amplihcation in a vicinity of corotation non-zero shear layer thickness is necessary 
(Savonije & Heemskerk 1990), the centrifngal instability does not develop, if we pnt < 
flex (Morozov 1977, Morozov 1979, Fridman 1990, Nezlin & Snezhkin 1990), the Kelvin- 
Helmholtz instability in this case is stabilized at M* = R{flex — flin)/cs < Merit = 8 ^/^ 
(Bazdenkov & Pogntse 1983, Torgashin 1986) (the value Merit = 8 ^/^ corresponds to flat 
vortex sheet stabilization, Landau 1944), and interaction of waves with opposite signs of 
energy at such M* values in absence of a central wall or density jump results to spontaneous 
radiation of neutral {Imu = 0 ) sound waves by jump fading at r —0 and r ^ 00 , and to 
considered instability if the wall is present. 

In view of told, the system of linearized equations of gas dynamics it is possible to 
reduce to system (Morozov 1977): 


dp 

dr 


2mfl 

ru 


p + Po{u‘^ 




( 1 ) 


^ = (^-^)Tb-(^ + -)e. (2) 

dr \Cg r^ J \ rut r J 

where u = u — mft is disturbances frequency with respect to Doppler shift, = 20(20 -|- 
rdft/dr) is squared epicyclic frequency, p{r) and ^(r) are peak functions of perturbed 
pressure and radial Lagrange displacement accordingly. The correspondence between the 
latter and radial component of perturbed speed is set by expression Vr = d^/dt — —iut^. 

In the area of uniformity the solution of system (1), (2) is a linear combination of 
modihed Bessel functions of complex argument of the integer order m: 


p{r) 


f^fmikinr^ B Krni^kinr^, Rp 7 T 7 Rqj 

CK^{k ), r > Rn, 


( 3 ) 


where 

kf = [40^ -{u- mfl‘^)]/cl. 


( 4 ) 
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Here the index i takes “in” and “ex” valnes. In (3) limitation of the solntion on inhnity 
is disconnted (for nniqneness of solntion choice we take Reki > 0 hereafter). 

On a velocity jnmp snrface bonndary conditions for the solntions (Morozov 1977) 
shonld be hold: 

^{r = Rci-0) = C{r = Rq + O), (5) 

p{r = R^-0)+ poQi^RciCiRci) = P{r = + 0) + pqVlI^Ru^{Rcl) , (6) 

and non-flow condition ^{Rp) = 0 on a wall snrface. 

Conpling the solntion for areas r < Rq and r > in view of conditions (5), (6) and 
condition at Rp and writing ont a non-flow condition of received system — the eqnality to 
zero of its determinant, — it is possible to receive the dispersion eqnation of considered 
model: 


^ex 

+ MlaS’ 

0 


(n) 
q; • ^ 


PI 


(n) 






- Min'S’ 

jip) 

^in 


-FPl 


= 0 , 


(7) 


where the following designations are introdnced: 


Rfl^i p _ Im{kinRp')RmikinRci) 

Cs Im(kinRil)PmikinRp') 


(j) _ , Fmi^jPj) 

* u — mVti * ^ Km{kiRj) ’ 

(j) _ 2mfli ^ xsR 

* u — mfli * ^ Im{kiRj) 

Here the prime stand for differentiation of Bessel fnnctions on argnment and top index 
shows that the given combination is written on radins of velocity or density jnmp. 

The eqnation (7) describes permitted eigenfreqnencies and the positive Imu presence 
(growth rate) means instability. 

2. The asymptotic colutions 

First of all it is necessary to note, that the eqnation (7) is satished identically at 
freqnencies ui = Vtinim — 2) or uj 2 = Oea;(m — 2). These roots describe gyroscopic modes 
of flnctnations, having in one of media only an azimnthal componenet of the wave vector 
kip — ^ (as kin = 0 ai u = ui and kex = 0 at u = U 2 )- This canse their nentral character 
— a nonzero flow of wave energy on radial coordinate in both media and, hence, ki ^ 0 is 
necessary for considered instability . 

In a limit of disappeaing small radins of a wall: D = {Rfi — Rp)/RQ —1, nsing asymp¬ 
totic decomposition Im{z) and Km^z) at z ^ 0 (Handbook of Mathematical Fnnctions... 
1964), it is possible to rednce eqnation (7) to a form: 




^ex ^in^Q. T ^inPin Pin ^ex^Q, T ^ex^ex — 0- 


( 10 ) 


The eqnation (10) accnrate to designations coincides with the dispersion eqnation from 
Morozov (1977), where a similar problem in absence of a wall was considered. 

In an incompressible limit {\kinRQ.\ ^ 1; \kexRQ.\ ^ 1) from (7) follows the solntion: 


U) ~ 


i^{Tlex ^inP T iS^ex ^in') T f ^ex \ \/(1 T Xfl) 11 -|- Tfl\ [ 


1 -/ 


( 11 ) 
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where I ~ [(1 — il>)^"^ + l]/[(l — —1], describing development of the Kelvin-Helmholtz 

instability. If not has places 1, (11) is rednced to asymptotic, received before (Moro¬ 
zov 1977): 




m{Qin + O ex ) + ex ^in - 1| 


( 12 ) 


With redaction D the presence of wall close to velocity jnmp snppresses this instability 
— growth rate monotononsly decreases and at ^ 0 takes place u mVtin i ■ 0. This 
resnlt, as well as (11), (12), does not depend on velocity jnmp “direction” on jnmp (i.e. it 
is valid both at Qin > ^ex and at Qin < ^ex)- 

In applied aspect a case of essential liqnid compressibility (|/cin-Rn| ni; \kexRii\ 
m) is mnch more interesting. Then (7) is rednced to a kind: 


kinRci 



kexRil 


th{kinRnD) ~ 0 . 


(13) 


When D 1 is not valid a solntion, received by Morozov (1977), valid only in a case 
Qin > flex and describing centrifngal instability, follows from (13): 



nr(Hin T flex) T 



( 14 ) 


However as against model considered by Morozov (1977), the presence at system of the 
aconstic screen (reflecting hrm wall) makes possible higher nnstable harmonics excitation 
(former at the flat jets stability analysis, snch harmonics were called “reflective” — see, 
for example, Payne & Cohn 1985) besides the basic un stable mode, described in different 
limits by expressions (11) and (14). Really, (13) admits the solntions of a kind: 


U) — OJq -|- 


( 15 ) 


where 

Cc’O 

dUJ — - 

DRU^ln - - 1 - iDRaiiVo - mnex)/cs ’ 


Wo - mQin y ’ 

n = 1,2,3,... is harmonic nnmber. The expression (15) is valid if Qin/\uj — mfiin\ 

1 ; Vtf>x/\^ ~ 'm^ex\ ^ 1 and |hw|/wo 1. As against (14) the expression (15) describes 
nnstable roots both at > Ogj, and at 0^^ < ^ex and in essentially snpersonic {R^\Qin — 
f2ex|/cs 3> 1) case, and is similar to analogons asymptotic of reflective harmonics in a 
plainparallel flow (see Morozov, Mnstsevaya & Mnstsevoj 1991). 

Finally in a strong compressibility limit again {\kinRn\ ^ m; \kexRn\ 3> m) bnt at a 
backlash between jnmp and wall so small, that \kinRnD\ 1, solntion describing another 
resonant mode follows from (13): 


U ~ mVtex + f 1 




SJL) 


Rq. - Rf. 


40L. 


(16) 


It is easy to note that first term nnder sqnare root sign, resnlting to instability in basic, 
represents the inverted r un time of a sonnd wave between jnmp and wall. Second term in a 
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square bracket defines stabilizing or destabilizing impact of difference of centrifugal force 
density on jump. Last term describes stabilizing influence of gyroscopic effects. 

Comparison (16) with results of Morozov & Mustsevoy (1991) and Morozov, Must- 
sevaya & Mustsevoj (1991) allows to conclude that (16) is long-wave (along a waveguide 
layer) limit for a harmonic with n = 1 described in a short-wave limit expression (15). 

3. Numerical solution results 

On Fig. 1.1-1.2 we bring dispersion curves, constructed on the data of the numerical 
equation (7) solution. For demonstration of considered ehect we have selected a case 
m = 2, Qin = 0, which is most “rigid” from the point of view of centrifugal stabilization of 
disturbances. With the purposes of facilitation of comparison with results of Morozov & 
Mustsevoy (1991) and Morozov, Mustsevaya & Mustsevoj (1991) the frequency normalized 
on sound z = uj{k^pCs) — ojRq,/ ijncs) is shown on figures. 

The basic results for a situation when gas between jump and wall is at rest are as 
follows: 

1) The presence of a central refiecting surface (acoustic screen) essentially infiuences 
on system stability. If in absence of this surface the vortex sheet is stabilized at M > 8^/^ 
(Bazdenkov & Pogutse 1983, Torgashin 1986), in a considered case the instability takes 
place even at M ^ 1. 

2) To change of parameters M and D for each azimuthal mode m there is the alter¬ 
nation of stability and instability areas. 

3) For disturbances with small azimuthal number m the basic unstable mode (the 
Kelvin-Helmholtz one) is stabilized at M = M* < 8^/^. M* grows as m and at m ^ 1 the 
basic mode ceases to be stabilized with Mach number growth, as it takes place in case of 
fiat vortex sheet in vicinity of a refiecting surface (see Morozov & Mustsevoy 1991). The 
reason of the last ehect as already was specified is instability of an acoustic resonance type. 

4) As well as predicts analytical equation (7) investigation, alongside with the basic 
mode the higher refiective unstable harmonics (see (15)) are excited and their quantity 
grows as m. 

5) Diherences from a case investigated by Morozov, Mustsevaya & Mustsevoj (1991) 
for higher harmonics at m 3> 1 when the curvature ehects are insignificant, are as small 
as well as for the basic mode. 

In a case flin > ^ex the results can be generalized as follows: 

1) In an incompressible limit i?n(h^m — ^ex) ^ Cg the disturbances dispersion law will 
be well coordinated with (11). 

2 ) At supersonic velocity diherence on jump both centrifugal and resonant modes 
develop. 

3) The centrifugal instability mode growth rate changes with change D extremely weak 
in a wide range of this parameter and is described by (14) successfully; only at H 0 the 
growth rate quickly decreases up to zero. 

4) Resonant modes growth rates are much less than centrifugal one and become com¬ 
parable to it only at H —0. 

4- Conelusions 

We shall move a total by formulating the basic conclusions about an opportunity of 
a central refiecting surface infiuence on results of experiments on rotating shallow water 
(Nezlin & Snezhkin 1990, Morozov et al. 1984, Morozov et al. 1985): 

1) The resonant instability of a considered type could not appreciably ahect spiral 
structure centrifugal instability modeling on shallow water because for parameters of the 
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unit “Spiral” its development time is significantly greater. Namely, speed jump and liquid¬ 
making pipe radii ratio such that D > 0.8 (Nezlin & Snezhkin 1990, Morozov et al. 1984, 
Morozov et al. 1985). 

2) In experiments on supersonic jump stabilization at speed difference growing up 
to R^i{Qin — ^ex) > 2y/2gH (Antipov et al. 1983) the resonant instability should result 
unstable regime renewal at further difference increasing, as against a situation considered 
by Bazdenkov & Pogutse (1983) and Torgashin (1986). Thus it would be curiously to repeat 
experiments of Antipov et al. (1983) in a wider range of relative speed of a liquid layers. 
At the same time even at smooth Mach number increasing the system will quickly pass 
through narrow domination areas of different harmonics n of various azimuthal modes m, 
which characteristic development time {Rn-Rp)/^, so any ordered wave structures 
observation will be hardly possible. 

The results obtained by us for a considered model situation cannot be applied directly 
to the analysis of stability of astrophysical objects and represent especially academic in¬ 
terest. However generalizing told it is necessary to specify necessity of further research of 
influence of a central density gradient (Zasov & Fridman 1987) on stability of galactic gas 
disks as the offered mechanism of disturbances exciting is rather strong. 
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